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ABSTRACT 

We report on an analysis of Fermi-LAT data from four year of observations 
of the nearby radio galaxy Centaurus A (Cen A). The increased photon statistics 
results in a detection of high-energy (> 100 MeV) 7-rays up to 50 GeV from the 
core of Cen A, with a detection significance of about 44(7. The average gamma-ray 
spectrum of the core reveals evidence for a possible deviation from a simple power- 
law. A likelihood analysis with a broken power-law model shows that the photon 
index becomes harder above Ej, ~ 4 GeV, changing from Ti = 2.74 ± 0.03 below 
to T 2 = 2.09 ± 0.20 above. This hardening could be caused by the contribution 
of an additional high-energy component beyond the common synchrotron-self 
Compton jet emission. A variability analysis of the light curve with 15-, 30-, and 
60-day bins does not provide evidence for variability for any of the components. 
Indications for a possible variability of the observed flux are found on 45-day 
time scale, but the statistics do not allow us to make a definite conclusion in this 
regards. We compare our results with the spectrum reported by H.E.S.S. in the 
TeV energy range and discuss possible origins for the hardening observed. 



-2- 



Subject headings: Galaxies: active - galaxies: individual (Cen A) 



Introduction 



The promi nent radio galaxy Centaurus A (NGC 5128), at a distance of ~ 3.8 Mpc 
(1' ~ 1.1 kpc) ([Harris et al. 20101 ) . is the closest active galax y to Earth. Often regarded 
as a prototype Fanaroff- Riley Class I iFanaroff fc Rilev 1974) rad i o source and as a mis - 
aligned BL Lac- type object at higher energies (jMorganti et al. 19921 : IChiaberge et al. 20011 ) . 
its proximity has made it on e of the best-studied extragalactic objects over a wide range of 
frequencies (e.g. Ilsrael 19981 . for review). Cen A is known for a complex and extended 
radio morphology, with two gia nt outer lobes ext ending over ~ 10° and oriented primar- 



ily in the north-south direction ( jFeain et al. 20 111 ). Optical images reveal the bright host 
galaxy bulge (~ 5' bulge radius) and the famous, warped dark lane of gas, dust and young 
stars (~ 12' in east-west extension) which obscures the inner part of the galaxy. Chan- 
dra X-ray observations show a one-sided, kpc-scale (up to ~ 4.5 kpc in projection) jet 
composed of several bright knots and diffuse emission, while high-resolution radio VLBI 
observations have also resolved jet and counter-jet feat ures on sub-parsec scales into dis- 
crete components (IKraft et al. 20021 : iMiiller et al. 201 ll ). These and related observations 
suggest that Cen A is a non-blazar source with its jet inclined at a rath er large viewing an- 
gle ~ 50° and characte ri zed by moderate (r adio) bulk flow speed ~ 0.5c flTingay et al. 1998 



Hardcastle et al. 2003 



Mii 



( jMeisenheimer et al. 2007 



ler et al. 20111). Its bolo metric luminosity of L ~ 10 42 erg/s 



van der Wolk et al. 20101 ). accompanied by indications for the 



lack of a dust torus, is thought to be p owered by gas accre t ion onto a supermassiv e black 
hole of mass M BH ^ (3 - 12) x 10 7 M Q jMarconi et al. 20061 : ICappellari et al. 2009h . 



At MeV energies, Cen A has been observed with both OSSE (0.05-4 MeV) and COMP- 
TEL (0.75- 30 MeV) onboard t he Compton Gamma-Ray Observatory (CGRO) in the period 
1991-1995 Jsteinle et al. 1998h . An agreement of the OSSE spectrum with the COMP- 
TEL one in the tra nsition region around 1 MeV, and correlated variability has been found 
(jSteinle et al. 19981 ) . At higher energies, a marginal (3a) detection of gamma-rays from the 
core of Cen A was reported with EGRET (0.1-1.0 GeV), bu t due its large angula r resolu- 
tion the association with the core remained rather uncertain (IHartman et al. 19991 ) . Unlike 
the initial variability (month-type?) seen at lower energies, th e flux detected by EGR ET 
appeared stable during the whole period of CGRO observation (ISreekumar et al. 19991 ). 



At VHE (> 100 GeV) energies, Cen A has also been detected (with a significance of 5a) 
by the H.E.S.S. array based on observations in 2004-2008. The results show an average VHE 
spectrum compatible with a power law of photon index T = 2.73 ± 0A5 sta t ± 0.2 syst and an 
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integral flux F(E > 250 GeV) = (1.56 ± 0.67 stat ) x 10~ 12 cm- 2 s- ^Aharonian et al. 20091 ). 
No evidence for variability has been found in the H.E.S.S. data set, but given the weak signal 
no certain conclusions can be drawn. 

Fermi-LAT has reported the detection of high energy (HE, > 100 MeV) gamma- 



rays from both th e core (i.e., within 0.1°) and the giant radio lobes ( lAbdo et al. 2010b 



Abdo et al. 2010af ). The analysis of 10 months of data revealed a point-like HE emission 
region coincident with the position of the radio core of Cen A, the emission being well de- 
scribed by a power-law function with a photon index fh 2.7, quite similar to the one in the 
VHE regime. Also, no variability (on 15 d and 30 d time scales) has been found. A simple 
extrapolation of the HE power-law spectrum to the VHE regime, however, fails to account for 
the TeV core flux as measured by H.E.S.S., which could indicate the need for an additional 
contribution towards the highest energies. The giant lobes were detecte d with a significanc e 



of 5<t and 8a for the northern and the southern structure, respectively ( lAbdo et al. 2010al ). 
The HE lobe spectrum could be described by a power-law function extending up to 2 or 3 
GeV with photon indices of T « 2.6. A recent analysis of a three times larger data set has 
confirmed the existence of these HE lobes, but also sh owed that the HE emission extends 



well beyond the WMAP radio image (jYang et al. 20121 ) . 



The apparent lack of significant variability features at GeV and TeV energies has so 
far precluded robust inferences as to the physical origin of the core emission in Cen A. 
Unfortunately, the resolutions of current gamma-ray instruments is not sufficient to localize 
the gamma-ray emitting region(s) either: The angular resolutions of both the H.E.S.S. array 
( ~ 0.1°) and Fermi-LAT (0.1° -1°, depending on energy) correspond to linear sizes of the 
gamma-ray emitting region(s) of about 5 kpc or larger. This ~ 5-kpc-region contains several 
potential gamma-ray emitting sites such as the central black hole, the sub-pc- or the kpc-scale 
jet etc. Based on the reported results, one thus cannot distinguish whether the gamma-rays 
observed from the core in Cen A originate in compact or extended regions. This motivated 
us to have a new look on the core emission based on four year of Fermi-LAT data. 



2. Fermi-LAT Data Analysis 
2.1. Data Extraction 

Fermi-LAT on board the Fermi satellite is a pair-con version telescope des igned to detect 



high-energy 7-rays in the energy range 20 MeV - 300 GeV (lAtwood et al. 20091 ). It constantly 



scans the entire sky every three hours and is always in survey mode. 

For the present analysis we use publicly available Fermi-LAT ~ 4 yr data from 4th 
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August 2008 to 1st October 2012 (MET 239557417-370742403). We use the Pass 7 data 
and analyze them using the Fermi Science Tools v9r27pl software package. The entire 
data set was filtered with gtselect and gtmktime tools and retained only events belonging 
to the class 2, as is recommended by the Fermi/LAT science team0. To reject atmospheric 
gamma-rays from the Earth's limb, events with zenith angle < 100 deg are selected. The 
standard binned maximum likelihood analysis is performed using events in the energy range 
0.1-100 GeV extracted from a 10° region centered on the location of Cen A, which is referred 
to as 'region of interest' (ROI). The fitting model includes diffuse emission components and 
gamma-ray sources within ROI which are not associated w ith Cen A (the model file is created 
based on Fermi-LAT second catalog (jNoland et al. 201 ll ). In the model file, the giant radio 
lobes were modeled using templates from WMAP-k band observation of the source which 
is extracted from NASA's SkyView. Although our previous results showed th at the HE 
extension and the radio lobe regions do not perfectly match (jYang et al. 20121 ). this does 
not affect the central parts of relevance here. The background was parameterized with the 
files gal _2yearp7v6_v0. fits and iso_p7v6source.txt and the normalizations of both components 
were allowed to vary freely during the spectral point fitting. 



2.2. Spectral Analysis 

Initially the continuum gamma-ray emission of the core of Cen A is modeled with a 
single power law. The normalization and power-law index are considered as free parameters 
then the binned likelihood analysis is performed. From a binned gtlike analysis, the best-fit 
power-law parameters for the core of Cen A are 

fdN\ , a f E \- 2 - 69± °- 03 

= 2-73 ± 0.12 x 10~ 9 — — . 1 

\dE J P V ; ViOO MeVy v ; 

This corresponds to an integral flux of 

F 7 = (1.61 ± 0.06) x 10~ 7 photon cm" V 1 , (2) 

with only statistical errors taken into account. The test statistic (defined as TS = 2 (log L - 
log L ), where L and L are the likelihoods when the source is included or not) is TS = 1978 
above 100 MeV, corresponding to a ~ 44 a detectio n significance. The results are consistent 



with the parameters found in (lAbdo et al. 2010bf ). namely photon index T = 2.67 ± 0.08 



1 http://fermi. gsfc.nasa.gov/ssc/data/analysis/documentation/ 
Cicerone/ Cicerone_Data_Exploration/Data_preparation. html 
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(between 200 MeV and 30 GeV) and integral flux (1.50 ±0.37) x 10~ 7 ph cm^s" 1 above 100 
MeV (model B). Figure [1] shows the spectrum of the core of Cen A obtained by separately 
running gtlike for 12 energy bands, where the dashed line shows the best-fit power-law 
function for the data given in Eq. ([[]). For the highest energy bin (56.2-100 GeV), an upper 
limit is shown. The spectrum shows a tendency for a deviation from a single power-law 
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Fig. 1. — Average high-energy gamma-ray (>100 MeV) spectrum of the core of Cen A 
(black points - t his work) as comp ared to the one based on the initial 10 month data set 
(blue squares - lAbdo et al. 2010a ). The dashed black line shows the power- law function 
determined from the gtlike. The blue and the red line show power-law fits to the energy 
bands below and above E b ~ 4 GeV, respectively. 

model with respect to the data above several GeV. Indeed, a \ 2 fit of the power-law model 
to the data gives a relatively poor fit with \ 2 = 39.7 for 9 degrees of freedom (dof), and 
its probability is P{x 2 ) < 2 x 10~ 5 . In order to investigate this in more detail, the core 
spectrum is modeled with a broken power-law model and gtlike tool is retried. The best-fit 
broken power-law parameters are 

-rv 2 



dEJ 



'1.19 ±0.08) x 10 



-13 



BP 



E 
E~ b 



(3) 
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and 

F 7 = (1.67 ± 0.06) x 10~ 7 photon cm" 2 s _1 , (4) 

with Ti = 2.74 ± 0.02 and T 2 = 2.12 ± 0.14 below and above E b = (4.00 ± 0.09) GeV, 
respectively. In order to compare the power-law and the broken-power-law model, a log 
likelihood ratio test between the models is applied. The test statistic is twice the difference 
in these log-likelihoods, which gives 9 for this case. Note that the probability distribution 
of the test statistic can be approximated by a x 2 distribution with 2 dof, corresponding to 
different degrees of freedom between the two functions. The results give P{x 2 ) — 0.011, 
which again indicates a deviation from a simple power-law function. The results of the data 
analysis with a broken power-law model reveal a hardening of the (average) gamma-ray core 
spectrum towards higher energies. The "unusual" break at 4 GeV could most naturally be 
explained by a superposition of different spectral components. In order to study this deeper, 
we divide the data set into two parts, i.e., (0.1- 4) GeV and (4-100) GeV. (Note that the 
4 GeV- value is obtained from binned maximum likelihood analyses). The core spectrum 
of Cen A in both energy ranges is then modeled with a power-law function and the gtlike 
tool is separately applied to these two energy bands. The photon index and flux between 
100 MeV and 4 GeV are I\ = 2.74 ± 0.02 and F 7 = (1.68 ± 0.04) x 10~ 7 photon cm-V 1 , 
respectively, and the test statistics gives TS=1944. The result is shown with a blue line in 
Fig. [TJ On the other hand, for the energy range (4-100) GeV we obtain T 2 = 2.09 ± 0.2 
and F 1 = (4.20 ± 0.64) x 10 -10 photon cm~ 2 s _1 , respectively, and a TS value of 124.4, 
corresponding to a w 11 o detection significance. This component is depicted with a red line 
in Fig. [TJ 



3. Temporal Variability 

Variability, if present, could provide important constraints on the emitting region(s). 
An observed HE flux variation on time scale t var , for example, would limit the (intrinsic) 
size of the gamma-ray production region to R' < -r^-ct var . However , previous HE and VHE 

\--\-Z i 



gamma-r ay observations of the core of Cen A with Fermi-LAT (lAbdo et al. 2010b! ) and 



H.E.S.S. ( lAharonian et al. 20091 ) did not find evidence for significant variability. Here we 
investigate whether the longer (4 yr) data set employed changes this situation. We thus 
divide the whole data set (from August 4th 2008 to October 1st 2012) into different time 
bins and generate light curves using the unbinned likelihood analysis with gtlike. Due to 
limited photon statistics the shortest time scale that one can probe is 15 days. In our 
analysis we generate light curves in 15, 30, 45 and 60 day bins. The normalization of the 
core and background point sources are treated as free parameters, but the photon indices of 
all sources and the normalization of the lobes are fixed to the values obtained in 100 MeV- 
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100 GeV energy range for the whole time period. Since no variability is expected for the 
underlying background diffuse emission, the normalization of both background components 
is fixed to the values obtained for the whole time period. To search for variability, a x 2 



2.5 - 




1.0 - 



54 800 55 000 55 200 55 400 55 600 55 800 56000 56200 

time (MJD) 

Fig. 2. — Gamma-ray light curve from August 4th 2008 to October 1st 2012. The bin size 
is 45 day. The background diffuse emission (both galactic and extragalactic) is fixed to the 
best-fit parameters obtained for the overall time fit. While some variability may be present, 
limited statistics do not yet allow to make definite conclusions. 

test was performed. The result for the light curve with 15 day bins is x 2 /d.o.f. = 1.22 and 
the probability is P{x 2 ) — 0.07. For the light curves with 30 day and 60 day bins we find 
X 2 /d.o.f. = 1.37 and x 2 /d.o.f. = 1.32, corresponding to P(x 2 ) = 0.04 and P(x 2 ) = 0.127, 
respectively. These results are consistent with no variability. Interestingly however, a similar 
test for the light curve with 45 day bins gives in x 2 1 'd.o.f. ~ 1.61 and P{x 2 ) — 0.015, 
indicating a possible variability on 45-day time scale. Unfortunately, because of limited 
statistics, we cannot make a definite conclusion in this regard. The light curve with 45 day 
bins is shown in Fig. |2l with the dot-dashed line indicating the flux from the source for the 
whole time period (result of likelihood analysis). 



4. Discussion and Conclusion 



In the case of high-frequency-peaked BL Lac objects, homogeneous leptonic synchrotron- 
self-Compton (SSC) jet models often provide reasonable descriptions of their overall spectral 
energy distributions (SEDs). For Cen A, however, classical one- zone SSC models (under the 
proviso of modest Doppler beaming) are unable to satisfactorily account for its core SEP up 
to the highest energies (cf. IChiaberge et al. 200ll ; iLenain et al. 20081 ; lAbdo et al. 2010bf ). It 
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seems thus well possib le, that an additional component contributes t o the observed emission 
at these energies (e.g.. iLenain et al. 2008r iRieger fc Aharonian 20091 ). The results presented 
here indeed provides support for such a consideration. Our analysis of the 4 yr-data set 
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Fig. 3. — Gamma-ray spectrum for the core of Cen A from high (Fermi-LAT , this work) 
to very high (H.E.S.S., blue squares) energies. The blue bowtie represents a power-law with 
photon index 2.74, and the red bowtie a power-law with photon index 2.09. The dashed lines 
show extrapolations of these models to higher energies. The power-law extrapolation of the 
low-energy component (blue lines) would under-predict the fluxes observed at TeV energies. 

reveals that the HE core spectrum of Cen A shows a "break" with photon index changing 
from ~ 2.7 to ~ 2.1 at an energy of Ef, ~ 4 GeV. This break is unusual in that the spectrum 
gets harder instead of softer, while typically the opposite occurs. For a distance of 3.8 Mpc, 
the detected photon flux F 1 = (1.68 ±0.04) x 10 -7 photon cm _2 s _1 for the component below 
4 GeV corresponds to an apparent (isotropic) 7-ray luminosity of L 7 (0.1 — 4 GeV) ~ 10 41 
erg s _1 . The component above 4 GeV, on the other hand, is characterized by an isotropic 
HE luminosity of L 7 (> 4GeV) ~ 1.4 x 10 40 erg s" 1 . This is an order of magnitude less when 
compared with the first component, but still l arger than the VHE lu minosity reported by 
H.E.S.S. L 7 (> 250 GeV) = 2.6 x 10 39 erg s" 1 fkharonian et al. 2009^ . All luminosities are 
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below the Eddington luminosity corresponding to the black hole mass in Cen A; nevertheless, 
they are still quite impressive when compared with the other nearby radio galaxy M87 
containing a much more massive black hole. 

Figure [3] shows the gamma-ray spectrum for the core of Cen A up to TeV energies. 
As one can see, the flux expected based on a power-law extrapolation of the low-energy 
component (below the break) clearly falls below the TeV flux reported by H.E.S.S.. Although 
the uncertainties in the photon index are large, it is clear that the spectrum becomes harder 
above 4 GeV. Remarkably, a simple extrapolation of the second (above the break) high- 
energy component to TeV energies could potentially allow one to match the average H.E.S.S. 
spectrum. These spectral considerations support the conclusion that we may actually be 
dealing with two (or perhaps even more) components contributing to the HE gamma-ray 
core spectrum of Cen A. Our analysis of the HE light curves provides some weak indication 
for a possible variability on 45 day time scale, but the statistics are not sufficient to draw 
clear inferences. 

The limited angular resolution (~ 5 kpc) and the lack of significant variability intro- 
duces substantial uncertainties as to the production site of the HE gamma-ray emission. 
In principle, the hard HE component could originate from both a very compact (sub-pc) 
and/or extended (multi-kpc) region(s). The double-peaked nuclear SED of Cen A has been 
reaso nably well-modeled up to a few GeV in term s of SSC processes occurring in its inner jet 
(e.g., IChiaberge et al. 2001t lAbdo et al. 2010bl ). In this context, the hardening on the HE 



spectrum above 4 GeV would indeed mark the appearance of a physically different compo- 
nent. This additional component could in principle be related to a number of different (not 
mutuall y exclusive) scenarios, suc h as (i) non-thermal processes in its black hole magneto- 
sphere (IRieger fc Aharonian 20091) . (ii) multiple SSC-emitting components (i.e., differential 
beaming) (iLenain et al. 20081) or (iii) ph oto-meson interactions of protons in the inner jet 
( Kachelriefi et al. 201o[ Sahu et al. 20121). Civ) 7-ray i nduced pair-cascades in a torus-like re- 
gion (at ~ 10 3 r s ) (e.g. iRoustazadeh fc Bottcher 2011 ) (v) secondary Compton up-scattering 
of host galaxy sta rlight (IStawarz et al. 200 6) or (vi) inverse-Compton (IC) processes in the 



kpc-scale jet (e.g. lHardcastle fc Croston 20111 ). What concerns the more compact scenarios 



(i)-(iv) just mentioned: Opacity considerations do not a priori exclude a near-BH-o rigin, but 



could potentially affect the spectrum towards highest energies (e.g. IRieger 20111 ). A SSC 
multi-blob VHE contribution, on the other hand, requires the soft gamma-rays to be due to 
synchrotron instead of IC processes, in which case correlated variability might be expected. 
Photo-meson (py) interactions with, e.g., UV or IR background photons (n 7 ) require the 
presence of UHECR protons, which seems feasible for Cen A. However, as the mean free 
paths A ~ l/(a pl n^K p ) of protons through the relevant photon fields are comparatively 
large, usually only a modest fraction of the proton energy can be converted into secondary 
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particles. Models of this type thus tend to need an injec tion power in high -energy protons 
exceeding the average jet power of ~ i0 43-44 erg/s (e.g. lYang et al. 20121 ). The efficiency 



of IC-supported pair cascades in Cen A, on the other hand, appears constrained by low 
accretion modes and the possible absence of a dust torus. Considering the more extended 
scenarios (v)-(vi): Partial absorption (~ 1%) of nuclear gamma-rays by starlight in the inner 
part of the host galaxy, and subsequent up-scattering of starlight photons could potentially 
introduce another HE contribution. However, the efficiency for this process is low, so that a 
high VHE injection power into the ambient medium is required, and the predicted spectral 
shape does not seem to match well. Compton-upscattering of starlight photon by energetic 
electrons in the kpc-scale jet also seems to have difficulties in reproducing the noted HE 
characteristics. 

Finally, let us mention that gamma-ray production may perhaps also be related to 
relativistic protons interacting with the ambient gas in the large (kpc) scale regions, e.g., the 
overall elliptical galaxy NGC 5128 or the densest part of its dust lane. Note that the 7-ray 
luminosity ~ 10 41 erg/s above 100 MeV is larger by two orders of magnitude than the 7-ray 
luminosity of the Milky Way, which could be related to a higher rate of cosmic-ray production 
and a more effective confinement in the case of NGC 5128. Moreover, gamma-rays might also 
be produced in a diluted -Rhaio ~ 30 kpc (halo) region of this galaxy. Despite the low density 
of gas, gamma-ray production on characteristic timescale t pp ~ 3 x 10 9 (n/10 _2 cm -3 ) -1 yr can 
be effective, even for a relatively fast diffusion of cosmic rays in this region. More specifically, 
the efficiency could be close to one, if the diffusion coefficient at multi-GeV energies does 
not exceed D ~ #LioApp ~ 1q29 cm 2 /s. This seems an interesting possibility, especially 
for the second (hard) HE component with photon index close to 2.1, in the context of its 
similarity to the gamma-ra y spectrum of the so-called 'Fermi Bubbles' around the center of 



our Galaxy ( [Su et al. 201)1 ). The much higher luminosity (by ~ 2 — 3 orders of magnitude) 
of the second component compared to the gamma-ray luminosity of the Fermi Bubbles seems 
quite natural, given the much larger energy available in Cen A, in particular in the form of 
kinetic energy of its jet. 

The results presented here provide observational evidence for an additional contribution 
at the highest energies and a more complex spectral gamma-ray behavior than previously 
anticipated. While considerations like those mentioned above may lead one to favor one 
production scenario over the other, none of them cannot be easily discarded. In fact, it is well 
conceivable that several of them contribute to the observed gamma-ray emission. Definite 
progress in this regard could be achieved in case of a significant detection of gamma-ray time 
variability. 
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